Defects in Graphene-Based Twisted Nanoribbons: 
Structural, Electronic and Optical Properties 
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We present some computational simulations of graphene-based nanoribbons with a number of 
half-twists varying from to 4 and two types of defects obtained by removing a single carbon 
atom from two different sites. Optimized geometries are found by using a mix of classical-quantum 
semiempirical computations. According with the simulations results, the local curvature of the 
nanoribbons increases at the defect sites, specially for a higher number of half-twists. The HOMO- 
LUMO energy gap of the nanostructures has significant variation when the number of half-twists 
increases for the defective nanoribbons. At the quantum semiempirical level, the first optically 
active transitions and oscillator strengths are calculated using the full configuration interaction (CI) 
framework, and the optical absorption in the UV/ Visible range (electronic transitions) and in the 
infrared (vibrational transitions) are achieved. Distinct nanoribbons show unique spectral signatures 
in the UV/ Visible range, with the first absorption peaks in wavelengths ranging from the orange 
to the violet. Strong absorption is observed in the ultraviolet region, although differences in their 
infrared spectra are hardly discernible. 



PACS numbers: 

I. INTRODUCTION 

Most recent advances in nanoscience and nanotechnol- 
ogy were made through the fabrication and characteri- 
zation of novel carbon-based nanostructures such as the 
buckminsterfuUerene , carbon nanotubes |2J [31 2] and 
graphenes [5J El |7]- Recently it is an indisputable fact 
that carbon has unique chemical properties among all el- 
ements of the periodic table, being essential to construct 
those nanomachines that existed in the natural world be- 
fore man appearance on Earth and work inside the cells 
of all living organisms [8 . 

FuUerenes and carbon nanotubes can be formed by 
rolling up a graphene sheet around a sphere and a cylin- 
der, respectively. Graphene shows many interesting fea- 
tures, such as extraordinary electronic transport proper- 
ties [3 Uni HU [12] related to its Dirac-like band structure. 
Indeed, a graphene sensor was able to detect the adsorp- 
tion of a single gas molecule thanks to the high sensitiv- 
ity of its electrical resistance to local changes in carrier 
concentration [13 . Another structure that can be made 
from a single sheet of graphene is a closed ribbon (one 
can think of a single wall carbon nanotube as a very wide 



closed carbon ribbon with a small radius). It is possible 
to conceive a closed carbon nanoribbon with one or more 
half- twists along its length. In particular, the nanoribbon 
with a single half-twist has the same features of the fa- 
mous Mobius strip, after the German mathematician Au- 
gust Ferdinand Mobius (1790-1868) who discovered it. A 
Mobius strip is a nonorientable surface in the Euclidean 
space 3?^, which means that a two-dimensional object 
transported around the surface can return to the point 
where it started looking like its image reflected in a mir- 
ror. General equilibrium equations for physical twisted 
strips were recently obtained by Starostin and Van Der 
Heijden [IJ. 

Molecules with twists and Mobius topology have al- 
ready been investigated in the literature, from both the 
theoretical and experimental viewpoints. More than 
forty years ago, Heilbronner [TS] deflned the concept of 
a Mobius aromaticity for cyclic molecules, predicting the 
stability of Mobius aromatic hydrocarbons with 4n p- 
electrons. Semiempirical and first principles calculations 
were used to investigate Mobius annulenes and twisted 
cyclacenes [161 [17]. In particular, the calculations for cy- 
clacenes with one, two and three half-twists [17] revealed 
a localisation of the twist over 2-4 benzene rings. First 



principles calculations using the Density Functional The- 
ory framework predicted the existence of a annulene with 
the topology of a Mobius strip with two half- twists [TH] • 



In 2003, Ajami et al. presented a paper reporting the 
synthesis of a neutral Mobius aromatic hydrocarbon |19j 
by combining two systems. The first one consisted in 
a flat aromatic ring with p-orbitals perpendicular to its 
plane, while the second system was a curved aromatic 
ring with pyramidalized bonds and p-orbitals parallel to 
the plane of the ring. Both systems were fused chemically 
through pericyclic reactions and five isomers were ob- 
tained, one of them with Mobius topology. Nevertheless, 
the several methodological errors pointed in the work of 
Ajami et al. with respect to the assignment of Mobius 
aromaticity to their newfound molecule [20j - a critique 
that deserved a reply [2T] - no one contests that a Mobius 
topology for a molecule was achieved in it. Other exam- 
ples of Mobius structures obtained experimentally are the 
crystals of NbSea synthesized by Tanda et al. and 
a protein with Mobius topology and insecticidal activ- 
ity found in some plants |23j . A recent account on the 
design of molecules with Mobius features is presented in 
Ref. 24J. 



In a previous paper [25] , we presented results of clas- 
sical force field, semiempirical and first principles simu- 
lations for twisted graphene nanoribbons with up to 7 
half- twists, investigating their structural stabilities, elec- 
tronic structure and optical properties. In this work we 
go further by studying the structural features of twisted 
graphene nanoribbons with defects created by the re- 
moval of a single carbon atom from different sites along 
the ribbon width. This kind of defect was chosen by 
us due to its size and aspect ratio as the removal of a 
single atom (single vacancy) is the most expected defect 
to be present in the nanoribbons. Other defects such 
as Stone- Wales, which are precursors of fractures in car- 
bon nanotubes [MJ , could also occur but they should be 
less probable and of higher formation energy. Among 
the methods for studying the effect of defects in carbon 
nanostructures one can cite the theory proposed by Wu 
et al. |27l based on the interatomic potential, incorpo- 
rating the effect of curvature and bending moment for 
curved surfaces, which allows for the determination of 
constitutive relations involving the stress, moment, strain 
and curvature as functions of the interatomic potential. 
Computational simulations were carried out using classi- 
cal dynamics and the semiempirical formalism. We show 
how the presence of defects affects the curvature and 
other structural characteristics of the nanoribbons, their 
electronic states and their infrared and UV-VIS spectra. 
We also perform a comparative analysis of these physi- 
cal properties as functions of the number of half-twists in 
the graphene nanoribbons, searching for molecular signa- 
tures that could be helpful to experimental testing. 
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FIG. 1: (a) Nanoribbon extracted from a graphene sheet, 
showing its dimensions L and W . (b) The A-defect created 
by removing a single atom of carbon from the nanoribbon 
central line (dotted line) of C-C dimers. (b) B-defect created 
by removing a carbon atom from a secondary line (dotted 
line). 



II. COMPUTATIONAL METHODOLOGY 

All twisted nanoribbons studied here were based in a 
single rectangular strip of graphene with armchair and 
zigzag edges, as in Ref. [2S]. The strip dimensions, L 
and W , are obtained by counting the number of C-C 
dimers parallel and perpendicular, respectively, to the 
armchair and zigzag sides. For the present investigation 
we have chosen L = 39 and = 7, as shown in Fig. 
1(a), which correspond to a strip length of 8.26 nm and 
width of 0.727 nm. This strip is larger than the one 
used in our previous work, L = 29, W = 5 [25] to allow 
for the evaluation of two different defect sites, and it is 
consistent in size and aspect ratio with some recent ex- 
perimental reports |2SJ [53] [50] . The nanostructure has 
280 carbon atoms and 86 hydrogen atoms are added to 
passivate their dangling bonds - 6 hydrogen atoms at the 
zigzag edges are removed when the strip ends are joined 
- so the chemical formulae of the flat rectangular strip 
and the closed nanoribbons are C280H86 and C28oH8o- 
A single defect can be obtained by removing one carbon 
atom from the strip. We have chosen to perform such 
removal from two sites denoted by the letters A and B, 
as depicted in Figure 1(b) and 1(c). The A-defect is ob- 
tained by deleting one carbon atom in the central line 
of C-C dimers parallel to the armchair border, while the 
B-defect is created by deletion of a single carbon atom in 
a secondary line. When the extremities of the strip are 
attached, the exact location of the defect along L does 
not affect the structural and electronic properties of the 
closed untwisted nanoribbons. For the twisted nanorib- 
bons the initial position of the defect site is random, the 
effects of arbitrary choice hopefully removed after suc- 
cessive geometry optimizations. 

Fifteen closed nanoribbons were built from the 
graphene strip, structures with the number of half-twists 
Nt — 0,1,2,3,4 and the number of defects Nd = 0,1, 
where the = 1 case takes into account the two dis- 
tinct defect sites A and B. Each nanoribbon was used 
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as input to classical molecular mechanics calculations in 
order to estimate the best structures to be used as start- 
ing geometries to compute in the quantum semiempirical 
framework. For this purpose we have selected the Forcite 
code and the COMPASS forcefield, which was optimized 
for condensed-phase atomistic studies, its parameters be- 
ing derived from ab initio calculations. To compute the 
molecular structures that minimize the classical total en- 
ergy we proceeded in three stages: first a geometry op- 
timization is made followed by a thermal annealing and 
new geometry optimization. In the first stage we start 
from crude structures obtained by simply twisting the 
graphene rectangle and joining its ends. These structures 
were optimized obeying some convergence thresholds: to- 
tal energy variation smaller than 2 x 10~^ kcal/mol, maxi- 
mum force per atom smaller than 10~^ kcal/mol/nm, and 
maximum displacement smaller than 10^^ nm. The op- 
timization algorithm uses a cascade of steepest descent, 
adjusted basis set Newton-Raphson and quasi-Newton 
methods. The non-bond interactions included electro- 
static and van der Waals terms, with an atom based 
summation method and cutoff distance of 3.2 nm. The 
optimized geometries correspond to local minima in the 
total energy hypersurface of each nanoribbon. The sec- 
ond stage was accomplished in order to achieve geome- 
tries corresponding to global minima (or improved lo- 
cal minima). In the thermal annealing calculation, each 
nanoribbon was submitted to a series of heating and cool- 
ing cycles to prevent structural trapping in a suboptimal 
energy minimum, a very likely outcome in the case of the 
optimized nanoribbons with a single defect and Nt > 0. 
Under such circunstances the defect site position along 
the length of the nanoribbon can be very important due 
to the smaller degree of symmetry in comparison with 
the Nt = case. Initial and mid-cycle temperatures 
were chosen to be 300 K and 1000 K, and 100 cycles 
were computed with 10 heating ramps per cycle and 100 
steps per ramp. After each cycle, a geometry optimiza- 
tion using the same criteria for the first stage was carried 
out. The NVT ensemble was adopted with a time step 
of 1.0 fs and a Nose thermostat with Q ratio of 1.0. In 
the third stage, we have carried out for all nanoribbons 
a new set of computations in order to reinforce our confi- 
dence that the obtained geometries are not the results of 
local mimima and/or used algorithms. We have chosen 
the same protocol of our previous work |25j using reactive 
molecular dynamics (Brenner-TersofF type) with different 
starting velocities and temperatures. This incorporates 
some stochastic aspects in the tests. The obtained partial 
structures were then reoptimized to check whether the 
obtained geometries were consistent with the expected 
miminum energy configurations. 

The best geometries computed in the classical mechan- 
ics formalism were used to begin a series of semiempirical 
calculations. The Austin Model 1 (AMI) Hamiltonian 
implemented in the Gaussian03 code was chosen to per- 
form a new geometry optimization within the restricted 
Hartree-Fock approximation. There are reports in the 
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FIG. 2: Nanoribbons with geometry optimized at the semiem- 
pirical level: (a) untwisted nanoribbon {Nt = 0), (b) Mobius 
nanoribbon {Nt — 1), (c) Nt = 2 nanoribbon, (d) Nt = 3 
nanoribbon, and (e) Nt = i nanoribbon. The bottom part 
of the Figure shows the defect sites for (f) A'^f = 1 A-defect, 
(g) Nt = 1 B-defect and (h) Nt — 2 A-defect nanoribbons. 
The rings with 9 and 5 carbon atoms at the defect sites are 
highhghted. 



literature [SH [31 [Ml 131] indicating that both AMI and 
first-principles Hartree-Fock calculations show compara- 
ble results for cages made of carbon or silicon supporting 
the applicability of the AMI semiempirical hamiltonian 
to investigate all-silicon and all-carbon clusters. With 
respect to the electronic structure, first principles DFT 
(Density Functional Theory) calculations tend to over- 
estimate the electronic correlation, underestimating the 
main energy gap (HOMO-LUMO) of molecular systems 
in comparison with experiment and semiempirical com- 
putations. Indeed, semiempirical CI calculations include 
some improvements on estimating the electronic corre- 
lation energy in comparison with AMI, predicting en- 
ergy gaps between the ground state and the first excited 
state intermediary between the DFT and simple AMI es- 
timates (those tend to follow the trend to predict energy 
gaps larger than experiment observed in first principles 
Hartree-Fock methods, which do not take into account 
electronic correlation effects). Using molecular mechan- 
ics, semiempirical and higher-level DFT-B3U^P meth- 
ods, Xu et al. [35j have found the ground-state isomers 
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of large fuUerenes (more than 100 atoms). Comparisons 
between the semiempirical and DFT-B3LYP results show 
that AMI, PM3, and MNDO semiempirical hamiltonians 
are notably less accurate for the prediction of the relative 
energies. On the structural viewpoint, however, bond 
lengths in carbon nanotubcs calculated within semiem- 
pirical approaches compare reasonably well with bond 
lengths found using different DFT exchange-correlation 
potentials [SB]. So we are confident that the structural 
and electronic structure features presented in this work 
can be considered, at least, as general trends for the fea- 
tures of real twisted carbon nanoribbons. 

Singlet ground states were calculated taking into ac- 
count the 1760 (1754) electrons of the defect-free (sin- 
gle defect) nanoribbons. Normal modes and the in- 
frared spectra were obtained afterwards and checked for 
negative frequencies. Finally, the oscillator strengths 
to estimate the UV-VIS spectra were found using the 
C.I. (Configuration Interaction) wavefunction approach 
taking into account the complete active space with con- 
tributions from 24 molecular orbitals (HOMO-11 up to 
LUMO+11). In the end, the C.I. ground state plus the 
first 127 excited states were acquired and analyzed. The 
final structures, after the semiempirical geometry opti- 
mization of the closed nanoribbons, are presented in Fig. 
2. 



III. STRUCTURAL PROPERTIES 

A few of the defect-free nanoribbons that minimize the 
semiempirical total energy exhibit some degree of symme- 
try. The Nt = nanoribbon has symmetry point group 
D2Qh- The Mobius nanoribbon {Nt = 1), on the other 
hand, has symmetry group Ci, differing from the corre- 
sponding atomic configuration used to start the calcula- 
tions, which belongs to the C2 point group. The opti- 
mized Nt — 2 nanoribbon has symmetry point group C2 
(starting from a C2 configuration), while for the Nt = 3 
(starting configuration with symmetry group C3) and 
Nt = 4: (initial symmetry C2) cases the symmetry group 
is C'l. The loss of symmetry in molecular strips is a 
consequence of the localization of the twists due to the 
search for a total energy minimum, being previously re- 
ported by Martm-Santamaria and Rzepa |17j and in our 
previous work [25 , both using the semiempirical frame- 
work. We believe that such twist localization is a general 
topological feature of all molecular Mobius systems. 

Nanoribbons with A and B defects exhibit structural 
changes in the carbon-carbon connectivity at the defect 
sites. The usual hexagonal rings of carbon atoms are re- 
placed there by two adjacent rings with different sizes, 
the biggest one having nine carbon atoms, while the 
smallest one has only five carbon atoms, as seen at the 
bottom part of Fig. 2 for the Nt — 1, A,B-defect (f,g), 
and the Nt = 2, A-defect (h) nanoribbons. In A-defect 
sites, the 9-carbon ring has C-C dimers along both arm- 
chair edges of the nanoribbons with A^t = 0, 1,4, but not 



for the A^t = 2, 3 cases. For the B-defect nanoribbons, 
the 9-carbon ring has always a C-C-C-C-C pentamer at 
the edge of the nanoribbon, except for the A^ = 4 case, 
where only a C~C dimer stays at the edge. The localiza- 
tion of the 5-carbon ring is variable, being either along 
the central line of the nanoribbon (A-defect, Nt — 0, 1,4, 
B-defect, Nt = 2, 3) with no C-C dimer along its borders, 
or with one C-C dimer along the nanoribbon border (A- 
defect, Nt = 2,3, B-defect, Nt = 0, 1,4). The formation 
of these 9, 5-carbon rings involve some degree of pyrami- 
dalization at the defect sites and, as we will see later, to 
a very pronounced increase in the local curvature of the 
nanoribbons in comparison with their defect-free coun- 
terparts. 



TABLE I: Heat of formation (in kcal/mol) calculated within 
the semiempirical formalism for each nanoribbon. 



Nt 


Defect-free A-Defect B-Defect 





1232.67 


1405.90 


1362.23 


1 


1298.16 


1458.67 


1413.09 


2 


1323.39 


1545.99 


1510.54 


3 


1446.19 


1594.70 


1534.67 


4 


1565.88 


1705.08 


1696.53 



Under the energetic viewpoint, the heat of forma- 
tion (AHf) of the twisted nanoribbons increases with 
the number of half-twists Nt, as one can see from Ta- 
ble 1. The defect-free nanoribbons have smaller values 
of AHf when compared with the respective values for 
the A- and B-defect nanoribbons. Starting at 1232.67 
kcal/mol {Nt = 0), AHf increases up to 1565.88 kcal/mol 
{Nt — 4:) in an almost linear fashion as a function of Nt- 
For nanoribbons with up to 7 half-twists, it has been 
shown that the total energy can be approximated as a 
quadratic function of Nt [15] . The A-defect nanoribbons 
for the same Nt have larger values in comparison with 
the B-defect ones, being in general larger by about 30- 
60 kcal/mol. One possible explanation for this difference 
is that it is more energetically expensive to remove a 
carbon atom from the central line of C-C dimers in a 
nanoribbon than from a secondary line. A carbon atom 
in a secondary line is closer to the edge of the nanorib- 
bon, a place where the energy cost of removing a single 
atom is lower. However, for Nt = 4, this difference de- 
creases to less than 10 kcal/mol. It may be possible that 
for nanoribbons with number of half-twists higher than 
4 the difference in the heat of formation between A- and 
B-defect nanoribbons is very small. 

To investigate the effect of the half-twists and defects 
on the curvature of the closed nanoribbons, we have de- 
fined the curvature angle 6, measured along the central 
line of carbon atoms shown in Fig. 1. This is the angle 
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FIG. 3: (a) The curvature angle is measured along the central 
line of atoms of the nanoribbon for three adjacent atoms. 
Curvature angle along the central line for all nanoribbons: 
(b) Nt = 0, (c) Nt = 1 (Mobius), (d) Nt = 2, (e) A^^ = 3, 
(f) Nt = 4. Each atom at the vertex of the angle is a site. 
Solid lines stand for the respective defect-free nanoribbon. 
Dashed lines stand for the A-defect nanoribbon and dotted 
lines correspond to the B-defect nanoribbon. 

between the straight Hnes connecting pairs of contigu- 
ous carbon atoms, as shown in Fig. 3(a). For the Nt = 
nanoribbon, as one can see from Fig. 3(b), the defect-free 
configuration has a practically constant value of 9 along 
its central line of atoms, about 9° (as expected from a 
regular polygon with 40 sides). The A-defect nanorib- 
bon shows a noticeable peak of 9, reaching about 27°, 
followed closely by the B-defect nanoribbon, with 22°. 
So, the presence of defects tend to create a kink in the 
nanoribbon. This can be understood in terms of mechan- 
ical rigidity: defects lower the mechanical resistance of 
the nanoribbons to shear stresses created by the process 
of curving the graphene sheet to form closed strips. This 
process seems to be reinforced when the closed strips are 
twisted, as shown in Fig. 3(c,d,e,f). 

For the Nt = 1 nanoribbon (Mobius case) with no de- 
fects, 9 oscillates smoothly between two almost symmet- 
rical regions with pronounced curvature, 9max ~ 15°, and 
two regions of low curvature, one with 9jnin ~ 0.17°, very 
flat, and the other with 9min ~ 5.8°. The insertion of an 
A-type defect produces a very pronounced increase of 9 
in one of the high curvature regions, with 9max reaching 
about 51°. The B-defect has not such a strong peak for 
the curvature angle, reaching about 25° at most. The sec- 



ond high curvature region remains practically preserved 
for both the defect-free and the defective nanoribbons. 
The largest curvature angle in the case of an A-defect for 
Nt — 1 occurs because the deletion of a central carbon 
atom weakens the mechanical resistance of the nanorib- 
bon by creating an empty space (a 9-carbon atoms ring) 
that crosses its entire width (see Fig. 2(f)). The B-defect, 
on the other hand, leaves an empty space only in one of 
the sides of the nanoribbon, with an extra hexagonal ring 
of carbon atoms reinforcing its structural integrity (see 
Fig. 2(g)). 

This explanation becomes more cogent by looking to 
what happens with 9 when the nanoribbon has two half- 
twists. The Nt — 2 nanostructure presents a 9-carbon 
ring contiguous to a secondary 5-(6-) carbon ring at its A- 
(B-) defect site, these rings being parallel to the nanorib- 
bon's width. As the 5-carbon ring has approximately the 
same mechanical resistance to shear of the 6-carbon one, 
it is expected that for both A- and B-defect nanoribbons 
the behaviour of 9 along the central line of carbon atoms 
should be very similar, as shown in Fig. 3(d). The defect- 
free nanoribbon has two very flat, symmetrical regions 
with 9min ~ 0.89°, and two high curvature regions with 
9max ~ 18°. The A- and B-defect nanoribbons have very 
close peak values of 9, 50° and 48° respectively, at one 
of the high curvature regions observed for the defect-free 
nanoribbon. The second region of high curvature looks 
pretty much the same for all structures, as occurred for 
Nt = 1 (and occur also for the Nt = 3, 4). 

The Nt = 3 nanoribbon has its 9 profile described 
in Fig. 3(e). For the nanoribbon with no defect three 
maxima of 9 can be seen, two of them nearly symmet- 
rical, with 9max ~ 27°, and a secondary maximum with 
^max ~ 11°- Three minima of 9 exist, two symmetri- 
cal flanking the secondary maxima with 9min ~ 0.50° 
and one minimum at the opposite side of the nanorib- 
bon, 9„iin ~ 1.7°. The A-defect nanoribbon 9 has a peak 
reaching 9„iax ~ 56° nearly at the same site of one of 
the principal maxima of the defect-free structure. For 
the B-defect, two peaks can be seen in the same region 
with 9max ~ 38°. The A-defect site forms a structure 
with two rings, one with 9 carbon atoms and the other 
with 6 carbon atoms, across the width of the nanorib- 
bon, while the B-defect site has one 9-carbon ring, two 
6-carbon rings and one 5-carbon ring in the same dispo- 
sition, enhancing its resistance to shear in comparison to 
the A-defect case, elucidating the highest 9max observed 
for the later. 

Finally, for Nt = 4: the mirror symmetry of 9 about the 
middle site observed in all previous cases {Nt = 1,2,3) 
for the free-defect nanoribbons disappears, giving rise to 
a distinct pattern of alternate maxima and minima (see 
Fig. 3(f)). From there one can see that 9 increases from 
2.7° to 14°, than decreases to 1.4°, rises again to 29° and 
tends to return to its initial value of 2.7°, as we move from 
site 1 to site 20. This sequence repeats itself from sites 
21 to 40. The nanoribbons with A- and B-defects show 
the same general behaviour for 9, reaching practically 
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from about 0.136 nm to 0.146 nm close to the defect site 
and is kept practically constant - about 0.140 nm - else- 
where. So the presence of defects produce a C-C bond 
length variation of about 4.3% at most in the central line 
of atoms. For the B-defect nanoribbons, the bond length 
variations are smaller due to the localization of the de- 
fect in a secondary line, parallel to the central line of C-C 
dimers. The bond length for B-defect nanoribbons tends 
to decrease near the defect, and its variation is smaller for 
nanoribbons with A^t = 2, 3, 4. In the right side of Fig. 4, 
the distance to the nanoribbon's centroid, dc, is plotted 
for all structures. For the Nt — case, dc is constant in 
the defect-free nanoribbon and oscillates slightly in the 
B-defect one. The A-defect nanoribbon dc oscillates more 
remarkably near the defect site. For Nf = 1,2,3,4, dc os- 
cillates with maxima in the regions of highest curvature 
angle (compare with Fig. 3). The amplitude of oscilla- 
tion is larger for greater values of Nt- The defects does 
not seem to affect remarkably the behavior of dc, except 
for a small shift in the localization of its maxima, as can 
be seen clearly for Nt = 3, 4. Finally, the length of the 
central line of atoms Lcl decreases from 8.49 nm to 8.38 
nm for the defect-free nanoribbon as Nt increases from 
to 4. For the A-defect nanoribbon, Lql starts at 8.48 
nm {Nt — 0), increases to 8.53 nm {Nt = 1), decreases 
to 8.38 nm {Nt — 2), 8.33 nm {Nt = 3) and increases 
again to 8.49 nm {Nt = 4). In the B-defect nanoribbon, 
Lcl decreases from 8.45 nm {Nt = 0) down to 8.38 nm 
{Nt — 3) and then increases to 8.44 nm {Nt = 4). 



FIG. 4: Dimer bond lengths (left) and distances between atom 
sites located at the central line of the nanoribbon and the 
nanoribbon centroid (right). The nanoribbons have (a) Nt = 
0, (b) Nt = 1 (Mobius), (c) Nt = 2, (d) Nt = 3, (e) N = 4. 
Solid lines, dashed lines and dotted lines indicate the defect- 
free, A-defect and B-defect nanoribbons, respectively. 



the same O^ax, 76° and 77°, respectively. These are the 
highest values for 6 observed for all nanoribbons and one 
can assume that this is an effect of the high degree of 
strain in a nanoribbon with such a big number of half- 
twists and the consequent coiling of the nanoribbon into 
three dimensions [37j. 

Besides the structural analysis using the curvature an- 
gle, we have also performed measurements of the bond 
lengths of all C-C dimers along the central line of atoms 
for each nanoribbon, as well as the distance between each 
carbon atom in the central line and the nanoribbon's cen- 
troid. The results of these measurements are presented in 
Fig. 4. One can see that the bond lengths for A-defect 
nanoribbons oscillate more strongly in a certain region 
for each Nt, which corresponds to the respective defect 
site. More pronounced peaks in bond length can be seen 
for Nt = 0,1,2,3, while for iVj = 4 there is a smaller, 
flat peak. The bond length along the nanoribbon ranges 



IV. ELECTRONIC STATES, UV-VIS AND 
INFRARED SPECTRA 



Subsequent to the geometry optimization at the 
semiempirical level, the electronic eigenstates were cal- 
culated for each nanoribbon. In Fig. 5 one can see 
the HOMO (Highest Occupied Molecular Orbital) and 
LUMO (Lowest Unoccupied Molecular Orbital) isosur- 
faces corresponding to a wavefunction amplitude of ±0.02 
for Nt = 0,1. The Nt = 2,3,4 HOMO and LUMO or- 
bitals were also obtained, but for the sake of brevity we 
have focused only in the Nt = 0,1 cases as they are rep- 
resentative of the general features observed for these or- 
bitals in all configurations investigated in this work. 

The HOMO states present phase alternation along the 
width of each nanoribbon for the Nt = 0, 1 defect-free 
and A-defect structures. For the B-defect nanoribbon, 
the pattern of alternating phases of the HOMO orbital 
seems to be rotated with respect to its armchair edge. 
LUMO orbitals, on the other hand, exhibit phase alter- 
nance along the length of all nanoribbons, except for the 
B-defect nanoribbon, which behaves in a similar fashion 
to the observed for the HOMO orbital. Both patterns of 
phase alternance associated to the defect-free geometries 
match the observations we have made in a previous pa- 
per |2S] considering L = 29, W = 5 twisted nanoribbons. 
For Nt = 0, defect-free structures, the molecular orbitals 
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HOMO 




LUMO 




FIG. 5; HOMO (top) and LUMO (bottom) orbitals for the 
Nt — 0, 1 nanoribbons. The first and third (second and 
fourth) rows correspond to the A^t = (Nt = 1) case. The 
first column refers to the defect-free nanoribbons, while the 
second and third columns show the molecular orbitals for the 
A- and B-defect nanoribbons, respectively. 



are delocalized. The twisting process, on the other hand, 
tends to locahze the molecular orbitals in a region of the 
molecular strip. When A-,B-defects are inserted, distinct 
effects on the HOMO (LUMO) localization appear. In 
the case of an A-defect nanoribbon with Nt = 0, 1, both 
HOMO and LUMO states confine electrons outside the 
region of the defect. On the other hand, for the B-defect 
case, electrons in the frontier orbitals show some degree 
of localization close to the defect sites. As the B-defect 
is localized near the edge of the nanoribbon, we presume 
that the distinct behaviour of the HOMO and LUMO or- 
bitals is due to the stronger interaction - in comparison 
to the A-defect structure - between dangling bond states 
at the defect site with electronic states along the edge of 
the nanoribbon. 

Eigenenergies calculated within the semiempirical 
method for the HOMO and LUMO electronic states are 
presented in Fig. 6, top-left. The defect-free nanoribbons 
show small variation of these energies as Nt is switched 
from to 4, with an average value of about -7.45 eV for 
the HOMO and -2.03 eV for the LUMO. The LUMO en- 
ergy decreases almost linearly when Nt increases from 
to 4, while the HOMO energy oscillates somewhat, in- 
creasing from iV( = to A^t = 1, decreasing from Nt = 1 
to Nt — 2, increasing from Nt — 2 to Nt = 3 and de- 
creasing again from Nt — 3 to Nt — 4. The A-defect 
nanoribbons have a more pronounced variation of the 
HOMO and LUMO energies as the number of half-twists 
changes. For the HOMO states, the energy starts at - 
7.46 eV {Nt = 0), decreases slightly to -7.48 eV {Nt = 1) 



and increases by almost 0.16 eV, reaching about -7.32 
eV {Nt = 2). Than it decreases to -7.44 eV {Nt = 3) 
and again to -7.49 eV {Nt = 4). The LUMO energy 
behaves in a reverse way with respect to Nt in compar- 
ison with the HOMO energy values. It decreases when 
the HOMO increases and increases otherwise: -2.01 eV 
{Nt = 0), -2.04 eV {Nt = 1), -2.30 eV {Nt = 2), -2.16 eV 
{Nt = 3), and -2.07 eV {Nt = 4). Finally, for the B-defect 
structure we have some oscillation of the HOMO, LUMO 
eigenenergies as the value of Nt is altered. These oscil- 
lations seem to mimic the oscillations observed for the 
A-defect nanoribbon in reverse. So when the HOMO en- 
ergy for the A-defect nanoribbon is augmented by adding 
a twist, the corresponding value for the B-defect nanorib- 
bon tends to get diminished. For the HOMO energies 
the values are -7.30 eV, -7.27 eV, -7.44 eV, -7.39 eV, and 
-7.48 eV, for Nt = 0,1,2,3,4, respectively. LUMO ener- 
gies are -2.24 eV, -2.25 eV, -2.12 eV, -2.15 eV, and -2.17 
eV for Nt — 0,1,2,3,4 in this order. Therefore for all 
cases, except for Nt — 2, the HOMO (LUMO) energies 
for the B-defect nanoribbon are larger (smaller) than the 
corresponding values for the A-defect nanoribbon. 

The HOMO-LUMO gap Eg is shown in Fig. 6, top- 
right. As the number of half- twists Nt changes for the 
defect-free nanoribbon. Eg exhibits a very small varia- 
tion, with minimum value of 5.41 eV for Nt = 1, 34 and 
maximum value of 5.43 eV {Nt = 0,2), a difference of 
only 20 meV. This could be compared with the data ob- 
tained for L = 29,W = 5 nanoribbons [25], with HOMO- 
LUMO gaps of 4.65 eV, 4.58 eV, 4.66 eV, 4.56 eV, 4.55 
eV for Nt — 0,1,2,3,4, respectively, the difference be- 
tween the minimum and maximum values being 110 meV. 
It seems that by increasing the size of the nanoribbon, 
the HOMO-LUMO gap tends to increase and becomes 
somewhat unsensitive to the variations in the number of 
half- twists. 

When a single defect is taken into account, a very 
pronounced variation is observed for Eg{Nt) as Nt is 
switched from to 4. For the A-defect nanoribbons. 
Eg starts at 5.45 eV {Nt = 0), than decreases a little 
bit to 5.44 eV {Nt = 1). For Nt ^ 2 a, sharp decrease 
occurs. Eg becoming 5.02 eV. For Nt = 3,4, the HOMO- 
LUMO gap energy increases to 5.28 eV and 5.42 eV, re- 
spectively. The Eg variation is 430 meV between the 
minimum {Nt = 2) to the maximum {Nt — 0). The B- 
defect nanoribbon also shows a large variation of Eg{Nt) 
in comparison with the defect-free case: 300 meV be- 
tween Eg{Nt = 2) = 5.32 eV and Eg{Nt = 1) = 5.02 
eV. Other values of Eg are 5.06 eV {Nt = 1), 5.24 eV 
{Nt = 3), and 5.31 eV {Nt = 4). One can note that the 
minimum value of Eg for the A-defect nanoribbon occurs 
for the same number of half-twists Nt = 2 the B-defect 
nanoribbon Eg exhibits its highest value. 

To investigate the excited states we have performed a 
full configuration interaction (CI) computation for each 
nanostructure. In order to build the CI determinants, 
the CI active space was formed using 24 molecular or- 
bitals from HOMO-11 to LUMO-l-11. The singlet ground 
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state is obtained by filling the 600 (598) lowest energy 
eigenstates with two electrons per level, leading to a total 
of 1200 (1196) electrons for the defect-free (A-,B-defect) 
nanoribbons. A total of 127 excited states were then cal- 
culated and the oscillator strengths for optical transitions 
between the singlet ground state and these excited states 
were computed to obtain UV/ Visible absorption spectra. 
The bottom of Fig. 6 shows the CI gaps obtained by tak- 
ing the energy difference between the first excited state 
Tl and the ground state GS (lines with open symbols) 
and by taking the energy difference between the first ex- 
cited state with transition oscillator strength different 
from zero Al and the ground state (first allowed optical 
transition represented by the lines with solid symbols). 
In all cases, the first excited state was a triplet, with opti- 
cal transition from the singlet ground state forbidden by 
the spin selection rule. The GS ^ Tl gaps for the defect- 
free nanoribbons (solid line, open squares) vary from 2.49 
eV {Nt — 1, minimal) to 2.63 eV {Nt — 0, maximum), a 
difference of 140 meV. The Nt = 2,3 nanoribbons have 
GS ^ Tl transition energies very close to the Nt — 
case, 2.60 eV and 2.61 eV, respectively. The iVt = 4 GS 

Tl transition energy is 2.52 eV, closer to the Nt = I 
structure. On the other hand, the optically allowed tran- 
sitions GS Al (solid line, solid squares) have higher 
values, but follow the same qualitative trend of the GS 

Tl energies, decreasing from 3.00 cV {Nt = 0, corre- 
sponding to violet color with wavelength 410 nm) to 2.67 
eV {Nt = 1, corresponding to blue color with wavelength 
460 nm), increasing to 2.76 eV {Nt = 2), increasing again 
to 2.80 eV {Nt = 3) and decreasing to 2.69 eV {Nt = 4). 
The overall variation in energy is 330 meV between the 
Nt = 1 and Nt = nanoribbons. 

The deletion of a single carbon atom to form an A- 
defect or a B-defect enhances the variation of the GS 
-> Tl and GS ^ Al CI gaps. Looking to the GS ^ 
Tl transition for the A-defect nanoribbon (dashed lines, 
open squares), one can see the GS — *■ Tl energy begin- 
ning at 2.65 eV {Nt = 0), decreasing to 2.55 eV {Nt = 1), 
decreasing again (and sharply) to 1.55 eV {Nt — 3), 
and then increasing to 1.75 eV {Nt = 4) and 2.49 eV 
{Nt = 4). For the optically active transitions GS — > Al 
the energy gaps are 2.84 eV {Nt = 0, 1, violet color with 
wavelength of 440 nm), 2.17 eV {Nt = 2, yellow color 
with wavelength of 570 nm), 2.26 eV {Nt = 3, green 
color, wavelength 550 nm), and 2.73 eV {Nt — 4, blue 
color, wavelength 450 nm). When a B-defect is formed, 
the CI gaps as functions of Nt show a behavior contrary 
to the observed for the A-defect case: when the CI en- 
ergy of the first increases as one half-twist is added, the 
corresponding CI energy for the later decreases, and vice- 
versa, with the exception when one switches from Nt = 3 
to Nt = 4. The B-defect nanoribbon GS ^ Tl energy 
gaps are 1.75 eV, 1.59 eV, 2.52 eV, 1.97 eV, and 2.16 eV, 
for Nt = 0,1,2,3,4 in respective order. For the GS 
Al optically allowed transitions, we find 2.23 eV {Nt = 0, 
yellow color, wavelength 560 nm), 2.08 eV {Nt — 1, or- 
ange color, wavelength 600 nm), 2.67 eV {Nt — 2, blue 
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FIG. 6: Top, left: HOMO and LUMO energy levels for the 
Nt = 0, 1, 2, 3, 4 nanoribbons. Top, right: HOMO-LUMO en- 
ergy gaps. For both plots, the solid lines indicate the defect- 
free case, while the dashed lines and dotted lines refer to the 
A- and B-defect nanoribbons, respectively. Bottom: energy 
gaps calculated within the semiempirical Full CI approach, 
with the curves with open squares, circles and triangles cor- 
responding to the difference of energy between the ground 
state (singlet) and the first excited state (in all cases, triplet). 
Curves with solid squares, circles and triangles correspond to 
the first allowed optical transition from the singlet ground 
state to a singlet excited state. Solid, dashed and dotted 
lines depict, in respective order, the defect-free and the A-, 
B-defect nanoribbons. 



color, wavelength 460 nm), 2.39 eV {Nt = 3, green color, 
wavelength 520 nm), and 2.71 eV {Nt = 4, blue color, 
wavelength 460 nm). So it seems to be possible to control 
the color in the visible spectrum a nanoribbon absorbs by 
changing the number of half-twists [25] . The addition of 
defects, on the other side, seems to increase the range of 
wavelengths in the visible spectrum accessible to optical 
absorption in comparison with the defect-free geometry. 

After obtaining the oscillator strengths for all tran- 
sitions in our CI calculations, we have plotted the 
UV/ Visible absorption spectra for all nanoribbons, as 
shown in Fig. 7. A set of Lorentzian curves with am- 
plitudes proportional to the oscillator strengths for each 
GS — > excited state (ES) transition and fixed widths 
of 0.038 eV was used in the plots. The absorption 
curves for the defect-free nanoribbons are shown in Fig. 
7(a),(b),(c),(d),(e) for Nt varying from to 4. When 
Nt = (Fig. 7(a)), three remarkable peaks appear at 
2.98 eV (violet color, wavelength of 420 nm), 3.73 eV 
and 4.11 eV (the first energy within the near ultraviolet 
range and the second within the ultraviolet B range). A 
very small peak also appears for an energy of 4.86 eV 
(corresponding to a photon in the ultraviolet C range). 
Nanoribbons with one or more half-twists show a richer 
absorption spectrum, with more peaks in comparison 
with the Nt = case. For A^t = 1 (Fig. 7(b)), the first 
significant peak appears at 2.88 eV, followed by peaks at 
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FIG. 7: UV/Visible spectra of the nanoribbons esti- 
mated by using the oscillator strengths for electronic tran- 
sitions between the singlet ground state and 127 ex- 
cited states. (a),(b),(c),(d).(c): defect- free nanoribbons 
with Nt = 0,1,2,3,4, respectively. (a'),(b'),(c'),(d'),(e'): 
Nt = 0,1,2,3,4, respectively, for the A-defect nanorib- 
bons. (a"),(b"),(c"),(d"),(e"): A^t = 0, 1, 2, 3, 4 for B-defect 
nanoribbons. 



3.68 eV, 3.83 eV and 4.10 eV. When iV* = 2 (Fig. 7(c)), 

the main absorption peaks appear at 2.94 eV, 3.88 eV 
(with a series of four smaller peaks between the two first 
ones), and 4.06 eV. In the Nt = 3 nanoribbon (Fig. 7(d)), 
one can sec that the single absorption peak of smallest 
energy observed in the Nt = 0,1,2 cases is replaced by 
two peaks, one small at 2.81 cV and the other one, big- 
ger than the first at 2.98 eV. A single small peak appears 
at 3.39 eV followed by a series of close absorption peaks 
between 3.84 cV and 4.17 cV. At last, for A^t = 4 (Fig. 
7(e)), a single absorption peak appears at 2.86 eV and a 
structure of close peaks with almost the same height is 
revealed in the 3.84 cV 4.17 eV energy range. Smaller 
absorption peaks also can be noted at 4.53 eV and 4.65 
eV. 

The A-dcfcct nanoribbons have their absorption spec- 
tra shown in Fig. 7 (a'), (b'), (c'), (d'), and (e'). For the 
Nt = configuration (Fig. 7(a')), five peaks are clearly 
visible, two near 3 cV, two near 3.8 eV, and one near 
4.1 cV. The last peak, after a closer look, exposes itself 
as a pair of nearby peaks. Therefore it seems that the 
addition of the defect splits the three sharp peaks at 2.98 
eV, 3.73 eV, and 4.11 eV observed in the Nt = defect- 
free nanoribbon. Looking now to the same nanoribbon 
Nt = with a single B-defect (Fig. 7(a")), the splitting 
of absorption peaks is apparent only for the peaks around 
3 eV, and a very small (but yet noticeable) absorption 
peak occurs at 2.23 eV, an energy equivalent to a green 
light photon with wavelength of about 560 nm. A simi- 
lar peak also appears at 2.10 eV ((Fig. 7(b")), Nt = 1, 
yellow-orange color) and 2.38 eV ((Fig. 7(d")), Nt = 3, 
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FIG. 8: Top: IR spectra for the Nt = (solid line) and Nt = 
1 (dotted line) defect-free nanoribbons. Bottom: difference 
of frequencies between corresponding normal modes of Nt = 
1 and Nt = defect-free nanoribbons (circles). The zero- 
frequency was shifted for better visualization 



green color), but it is absent in the Nt = 2,4 B-defect 
nanoribbons ((Fig. 7(c"),(e")). In the A-dcfcct nanorib- 
bon with Nt = 2 (Fig. 7(c')), a very small absorption 
peak occurs for an energy of 2.18 eV (yellow color) and 
in the A^t = 3 case (Fig. 7(d')), we have another small 
peak at 2.25 eV. A common feature of all nanoribbons 
with defects is a larger number of peaks in comparison 
with their defect-free counterparts, but with smaller ab- 
sorption intensities. 

Next to the calculation of the electronic absorption 
spectra in the UV/Visible range, we now comment the 
infrared spectra calculated for the defect-free nanorib- 
bons with Ni = 0. 1 (Fig. 8). We do not show the 
infrared spectra for the other nanoribbons (defect-free 
with more half-twists, A-, B-defect) because they resem- 
ble very closely the spectra we arc about to discuss. In- 
deed, as shown in the top of Fig. 8, the infrared spectra 
for the Nt = and Nt = 1 nanoribbons are so similar that 
their differences can be seen only when we plot it using 
a logarithmic scale. The A^t = 1 IR peaks (dotted fine) 
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appear practically at the same energies of the Nt = 
case (solid line), with some low intensity peaks exclusive 
of the Nt — 1 nanoribbon observable only for the 200- 
400 cm-\ 590-640 cm'^, and 1100-1200 cm'^ frequency 
ranges. A comparison of the normal mode frequencies 
uj for each nanoribbon was carried out by plotting the 
difference between the corresponding normal mode fre- 
quencies uj{Nt = 1) — uj{Nt = 0) versus ui{Nt — 0), 
as shown in the bottom part of Fig. 8. A correspond- 
ing normal mode is defined as follows: first we enumer- 
ate in crescent order of frequency the normal modes of 
each nanoribbon; second, as both nanoribbons have the 
same number of vibrational eigenstates (1080 total), wc 
can make an one-to-one correspondence between them. 
For frequencies between and 1000 cm~^, the TVj = 1 
nanoribbon has normal modes with frequencies a little 
bit higher than the frequencies for the untwisted geome- 
try, with maximum difference of about 21 cm~^. In the 
frequency range 1000-1600 cm~^ this trend is reversed, 
with the Nt — I frequencies being in general lower than 
the Nt — ones, with minimum difference of about -27 
cm~^. Between 1600 and 1800 cm~^ and for the high- 
est frequency modes around 3100-3200 cm~^, the corre- 
sponding normal modes of both nanoribbons, Nt = 0,1, 
have very close frequencies. The mean value of the dif- 
ference uj{Nt = 1) — uj{Nt — 0) was -0.60 cm^^, with 6 
cm^^ of standard deviation. So the normal modes of the 
Nt — I nanoribbon tend to have slightly smaller frequen- 
cies than the corresponding normal modes of the Nt = 
structure. 



V. CONCLUSIONS 

Recent advances in the synthesis of molecules with 
half-twists and Mobius topology stimulate the investiga- 
tion of their structural, electronic and optical properties. 
At the same time, the astounding properties of graphene 
are driving huge research efforts worldwide. One of the 
possible nanostructures derived from graphene is pre- 
cisely a closed-loop twisted nanoribbon. In this paper 
we pointed out some results of classical and quantum 
semiempirical simulations of graphene-derived twisted 
nanoribbons with length L = 39 and width W — 7 (pa- 
rameters defined as in Ref. [25]) with defects created by 
the deletion of a single carbon atom in two distinct sites, 
comparing their characteristics with the corresponding 
defect-free nanoribbons. Hydrogen atoms were added to 
passivate dangling bonds (except at the defect sites). 

After obtaining the best geometries for the nanorib- 
bons using a combination of classical dynamics and clas- 
sical annealing, a set of closed nanoribbons with number 
of half-twists (Nt) varying from to 4, defect-free and 



with two types of defect (A and B) was submitted to 
a process of geometry optimization using the quantum 
semiempirical Hamiltonian AMI. In defective nanorib- 
bons, the pattern of hexagonal carbon rings is replaced 
at the defect location by two adjacent rings with 9 and 
5 carbon atoms (9-5 rings). The spatial disposition of 
the 9-5 rings in the nanoribbons seems to depend on the 
number of half-twists as well as on the type of defect and 
contribute, through the pyramidahzation of C-C bonds, 
to increase local curvature. Structural features of the 
carbon nanoribbons were evaluated and it was demon- 
strated that the insertion of defective sites significantly 
increases the local curvature angle 9. In general, A-defect 
nanoribbons have more local curvature at the defect site 
in comparison with the B-defect ones. C-C bond lengths, 
however, do not change meaningfully when one compares 
the defective and defect-free structures among themselves 
(less than 4% of difference). 

Frontier molecular orbitals are delocalized for Nt = 0, 
defect-free nanoribbons, and localized in the twisted 
(specially in the defective) nanoribbons. By inserting A- 
(B-) defects, the HOMO and LUMO states have larger 
amplitudes in regions apart from (at) the defect-site. The 
HOMO-LUMO gap of the A-defect nanoribbons increases 
(decreases) when the HOMO-LUMO gap of the B-defect 
nanoribbons decreases (increases) as functions of Nt . The 
same happens for the CI gaps. In the defect-free nanos- 
tructures, the first optically active transitions (with non- 
zero oscillator strength) involve photon energies corre- 
sponding to violet {Nt = 0) and blue {Nt = 1,2,3,4) 
colors. Defective nanoribbons, on the other hand, have 
the first optically active transitions with energies rang- 
ing from 2.08 eV (B-defect, Nt = 1, orange color) to 2.84 
cV (A-defect, Nt = 0,1, violet color). The most intense 
UV/ Visible absorption peaks were observed in the ultra- 
violet range, with smaller peaks in the visible spectrum. 
Contrasting the defect-free nanoribbons absorption lines 
with the A-, B-defect ones reveals a richer structure of 
peaks for the later. Adding half-twists also contributes 
to increase the number of peaks, in some cases through 
splitting (for example, the Nt = S defect-free nanorib- 
bon). Distinct nanoribbons have very similar infrared 
spectra, even in the presence of A or B defects. 
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